TT-Junction behavior and Andreev bound states 
in Kondo quantum dots with superconducting leads 
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We investigate tiie temperature- and coupling-dependent transport through Kondo dot contacts 
with symmetric superconducting s-wave leads. For finite temperature T we use a superconducting 
extension of a selfconsistent auxiliary boson scheme, termed SNCA, while at T = a perturbative 
renormalization group treatment is applied. The finite-temperature phase diagram for the O-tt 
transition of the Josephson current in the junction is established and related to the phase- dependent 
position of the subgap Kondo resonance with respect to the Fermi energy. The conductance of 
the contact is evaluated in the zero-bias limit. It approaches zero in the low-temperature regime, 
however, at finite T its characteristics are changed through the coupling- and temperature-dependent 
O-TT transition. 

PACS numbers: 74.50.-|-r, 72.15.Qm 



I. INTRODUCTION 

The physics of charge and spin transport through 
Kondo quantum dots is paradigmatic for interface prob- 
lems with strong correlations mediated by the contact. 
Novel effects, taking place in junctions with quantum 
dots between normal metal leads, have been intensively 
studied for a long periodiSiSi^ The Josephson current 
through a localized spin state was first considered by 
Shiba and Soda*^ Later on Glazman and Matveev in- 
vestigated more thoroughly the supercurrent through a 
single resonant state, as well as through a distribution of 
such impurity states^. Eventually, a Kondo quantum dot, 
which is coupled to a normal and to a superconducting 
lead, is a further notable system which allows to approach 
the interplay between Kondo effect and Andreev reflec- 
tionsii^ The progress in the miniaturization of electronic 
devices now makes the investigation of electronic trans- 
port through a single Kondo impurity technically feasi- 
ble. To date, several groups have reported on transport 
measurements of such nanoscale devicesi24ii*i2iiiiiiiSii&. 
Besides possible applications as, for example, the study of 
nonlocal spin-entangled pairs^, these quantum dot con- 
tacts are fascinating on fundamental grounds, because 
they are the most elementary realization of a "strongly 
correlated contact" . 

From the theoretical side it has been well apprehended 
that a phase-sensitive subgap state is formed, which is to 
be interpreted as a Kondo resonance, if the Kondo scale 
Tk is larger than the gap A of the superconducting leads, 
TkI ^ ^ 1- Andreev scattering processes induce a de- 
pendence of the subgap-state energy on the phase differ- 
ence of the superconductors in the leads. The inter- 
ference of Andreev scattering with Kondo-type spinflip 
processes leads to a non-trivial behavior of the Joseph- 
son current -phase relation /s((^) in Kondo quantum dot 
junctions. A transition from a 0-junction to a 7r-junction 



arises on account of the distinct nature of the spin ground 
states in the strong and in the weak coupling regime. In 
the strong coupling limit, T/^/A ^ 1, the ground state is 
a spin singlet due to the Kondo screening of the impurity 
spin, and the Coulomb blockade is lifted by the formation 
of the Kondo resonance. In this case, coherent Cooper 
pair transmission occurs without affecting the spin of the 
electrons in the pair (0-junction). In the opposite limit of 
weak coupling, Tjf/A <C 1, the Kondo screening is sup- 
pressed for temperatures T well below the critical tem- 
perature Tc of the superconductors as the Cooper pairs 
in the bulk cannot be broken for s-wave pairing symme- 
try. Then the ground state is a Kramers degenerate spin 
doublet, and a single subgap resonance of width ~ Tk 
in the impurity spectrum is formed, split off from the 
continuum spectrum. In this regime, retarded, coher- 
ent pair transmission is still possible, but for energetic 
reasons (no double occupancy of the dot) the temporal 
sequence of the transmitted electrons with opposite spin 
is reversed, leading to a 7r-shift in the current-phase re- 
lation (7r-junction)jiS see Ref . l2dl for a more detailed dis- 
cussion. Since for weak coupling the subgap state forms 
below and for strong coupling it moves above the Fermi 
energy!i2i2i the current-phase relation /s(0) is also re- 
lated to the position of the resonance. This behavior 
of a Kondo quantum dot should be contrasted to the 
case when the impurity state is not a dynamical quan- 
tity and its magnetic moment is fixed. The latter case 
is analogous to junctions with ferromagnetic interlayers, 
where Andreev subgap states are generated both below 
and above the Fermi energy, being split with respect to 
their spin polarizationi22iSi2i2^ 

Whereas the supercurrent through a Kondo correlated 
junction has been investigated successfully within several 
approaches, the conductance of the contact is much more 
difficult to study, as it involves the quasiparticle current. 
It is essential to distinguish between Kondo 'point con- 



2 



tacts and Kondo quantum dot devices. A Kondo impu- 
rity in a point contact or orifice introduces an additional 
scattering channel and tends to reduce the transmission 
similar to the Kondo effect in bulk metals. In contrast, 
for a Kondo quantum dot device, the quantum dot pro- 
vides the only transmission channel. For temperatures 
above Tk and a quantum dot energy level well below the 
Fermi energy, the transmission channel is "almost closed" 
as the charge tunneling is suppressed by the Coulomb 
blockade and because the level is off resonance. Through 
the formation of the Kondo resonance at temperatures 
below Tk, on- resonance tunneling enhances the trans- 
mission up to the quantum limit While this inverse 
relation between the two junction types is rather obvious 
for normal conducting leads, it has more profound con- 
sequences in the case of superconductors. In the present 
paper we study the intrinsic conductance which charac- 
terizes the quasiparticle current through Kondo dots be- 
tween two superconductors. The quasiparticle current is 
not related in a simple way to the supercurrent and the 
question arises, if the O-tt transition may already mani- 
fest itself in the zero bias conductance of the Kondo dot 
contact. Choi et alm^ also investigate a Kondo quantum 
dot with superconducting leads to calculate the Joseph- 
son current at T = 0. However, to determine the low-T 
conductance through the O-tt transition, they consider a 
Kondo quantum dot with an additional resistive shunting 
(resistively shunted superconducting junction, RSCJ) in 
the overdamped regime and compare the crossover for 
Tk ~ 0.5A with the measured conductance of gated car- 
bon nanotube quantum dots coupled to superconduct- 
ing Au/Al leadsS^. Although the RSCJ modelling may 
well apply to the considered experiments, it actually does 
not refer to a Kondo quantum dot as defined above. It 
does not analyze the quasiparticle current through the 
Kondo impurity but rather the phase slips of the super- 
current. The respective conductance Gs in the RSCJ 
model grows exponentially with the inverse temperature 
Gs/Gm ~ exp{hIs/eT) where Gjv is the conductance in 
the normal state. 

In this article we address two related subjects which 
are relevant for Kondo dots between two superconduct- 
ing leads. On the one hand, a finite-temperature phase 
diagram for the 7r-junction behavior has not yet been 
presented. It allows to identify the coupling strength at 
which the O-tt transition sets in, but it also renders the 
regime where the transition may be observed by vary- 
ing the temperature at fixed coupling strength. On the 
other hand, we calculate the intrinsic conductance of the 
Kondo dot contact — a quantity which is not related 
to the phase difference as directly as the supercurrent. 
Nevertheless, as both quantities, the supercurrent and 
the quasiparticle current, are controlled by the position 
of the subgap resonance, they both display a transition in 
their dependence on temperature and coupling strength, 
as we will discuss in this paper. For this investigation 
we will always consider the generic case of a symmetric 
junction and s-wave pairing symmetry in the leads. 



In Sec. ni^e briefly introduce the technique used for 
the Kondo correlated junctions, an extension of the non- 
crossing approximation (NCA) to the superconducting 
state (SNCA) . Some explicit details about the derivation 
of the SNCA, its evaluation and regime of validity, as 
well as the calculation of the supercurrent and the con- 
ductance are deferred to the appendices. Sec. IIIII expands 
the discussion of the O-tt transition in certain aspects be- 
yond what has been presented in the literature on this 
topic. Specifically, we focus on the temperature depen- 
dence of the transition and introduce a phase diagram. 
In Sec. IIVI we address the temperature- and coupling- 
dependence of the intrinsic conductance of the Kondo 
dot contact. 



II. KONDO IMPURITY BETWEEN TWO 
SUPERCONDUCTORS 

The system of conduction electrons in the left and right 
lead interacting with a single-channel magnetic impurity 
or quantum dot is modeled by an infinite-?/ Anderson 
Hamiltonian. The s-wave superconducting state in the 
reservoirs is treated within standard BCS mean-field the- 
ory. The complete Hamiltonian then takes the form (see 
Fig. ^for a graphical layout of the contact) 

H = Ho + Hbcs + Hqd + XQ (1) 

with 

kcr 

HbCS = - E (ckTaclkia + h.C.) , 
ka 

(T k(T 

Here we have adopted a slave-boson representation^^ for 
the dot states, where the local creation and annihilation 
operators for an electron in the dot (d-) orbital with 
spin a and energy Cd are decomposed as, for example, 
= /J6. The operators and create a singly oc- 
cupied or an empty occupied impurity state, whenever 
an electron hops onto or off the dot, respectively, and 
obey the canonical fermion and boson commutation re- 
lations. Their dynamics are restricted to the physical 
Hilbert space by the operator constraint 

cr 

which will be enforced exactly by taking the limit of the 
parameter X ^ oo (see Appendix A)i^ Moreover, c^^^ 
creates a conduction electron in the left (L) or right (R) 
superconductor, a = L, R. The hybridization of these 
electronic states in the leads with the quantum dot state 
is parameterized by Va- For convenience, we introduce 
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FIG. 1; Quantum dot coupled to two superconductors. Tl 
and r_B denote the effective couplings to the left and right 
lead, 4>L and (f>R label the phases of the corresponding super- 
conducting order parameter. The BCS gap A is assummed 
to be equal in both superconductors. In the model, defined 
by Eq. (0, the local Coulomb repulsion U is set to infinity. 



the effective couplings Fq = nNoV^ and F = F^ + F^, 
where Nq is the density of states at the Fermi energy in 
the normal conducting state. 

The BCS part of the Hamiltonian can be easily solved. 
The normal and the anomalous local advanced/retarded 
conduction electron Green's functions are defined 
as G^/^W = ±Ek*^(Tt)({Ck,,W,cLa(0)}) and 

= ±Ek*^(TO({ckTaW,c_kia(0)}>, respec- 
tively. The gap equations defining the order parame- 
ter Aa in the two superconductors are given by Aq = 

^Bcs^'ki'^-'kia'^kla) ■ Subsequent consideration 

the amplitude of is assumed to be equal on both sides, 
i.e. 



A, = |A| 



(3) 



For the local conduction electron density of states per 
spin and the corresponding anomalous contribution one 
obtains. 



(4a) 



..(^)^^^^^^^-Re4|^, (4b) 



where both spectral functions have been normalized to 
No. 

For the greater part of this paper the Kondo dynam- 
ics of the quantum dot at finite temperatures will be 
described within a selfconsistent approach, where the lo- 
cal gauge symmetry on the dot is preserved by means 
of conserving approximations, derived from a Luttinger- 
Ward functional'^ ^. We will use a generalization of 
the well-known non-crossing approximation (NCA)^2i2^ 
for superconducting leads, the "superconducting NCA" 
(SNCA), to include retarded Cooper pair tunneling. The 
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FIG. 2: Generating functional for the extension of the NCA 
to the broken-symmetry state (SNCA). The solid, wavy and 
dashed lines represent the conduction electron, slave boson 
and pseudofermion propagators, respectively. 



Luttinger-Ward generating functional $ for the SNCA is 
depicted in Fig. [3 The leading term of 0(T) in a selfcon- 
sistent expansion corresponds to the NCA (first diagram 
in Fig. |2J|. For normal conducting leads, the NCA is 
known to give a satisfactory, quantitative description of 
the spectral features in the case of infinite £/ j?"^i?^i?^i?7i?? 
in the absence of magnetic fieldf^^iiS and for tempera- 
tures down to T « 0.1 Tffi^. However, in the case 
of superconducting leads the NCA completely neglects 
Andrcev scattering contributions (Cooper pair tunnel- 
ing through the dot), which is crucial for the Joseph- 
son current and which will also induce significant renor- 
malizations of the normal quasiparticle current, as seen 
below. Therefore, the NCA is extended to include the 
next-to- leading term of order O(r^) (second diagram in 
Fig. 12), which contains two anomalous lead Green's func- 
tions, constituting the SNCA^i Similar, but simplified 
methods, employing an elastic scattering approximation, 
have also been used by Bickers and Zwicknag!^^ and by 
Borkowski and Hirschfeld^. A detailed discussion of the 
SNCA is deferred to Appendix A. It will be seen that 
the SNCA describes, to leading selfconsistent order, the 
coherent transmission of Cooper pairs via the formation 
of retarded Cooper pairs on the dot even though the exis- 
tence of equal-time Cooper pairs on the dot is prohibited 
by the local Coulomb repulsion. Superconducting Kondo 
dot junctions have recently been considered also within 
a mean field approach to the dot dynamicS(^^i^ which 
tends to overestimate the Cooper pair correlations on the 
dot due to the assumption of static rather than retarded 
pairs on the dot. 

From the generating functional in Fig. [21 a set of cou- 
pled integral equations for the pscudoparticle selfenergies 
can be derived as well as an expression for the local spec- 
tral function of the quantum dot and its corresponding 
anomalous part. These equations are solved numerically. 
The explicit expressions are discussed in Appendix fXl 

To calculate the Josephson current we use the formula 
first presented by Clerk and Ambegaokar— 

2e F f 
Is{<l>) = ———smcj) / dujf{uj) 



h nNo 



(5) 



which is rederived in Appendix ^ The quantities in this 
current relation are defined as follows: 4> = 4>l — (t>R de- 
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notes the phase difference between left and right lead, 
/(w) is the Fermi function and, for convenience, we 
extracted the explicit phase dependence from the off- 
diagonal Green's functions, 

•^<fM = cos(^0^fH 

where !F^\u)) and F^{u!) are the anomalous parts of 
the Green's function of the impurity d-level and of the 
conduction-electron Green's function in lead a, respec- 
tively (cf. Appendix IB|| . 

The zero bias conductance G = dI/dV\v=o is calcu- 
lated from the quasiparticle current in the limit of small 
bias^, 

G = -2^r J diu pM \mg^{u:) , (6) 

with gf{uj) the normal part of the impurity Green's func- 
tion. As we consider a symmetric coupling to the two 
leads with equal spectral densities (Fl —Tji = ttNqVI ^, 
F EE Fi -I- F_R = 2Fl,_r, and = Pr{uj) = p{uj)), all 

contributions with anomalous as well as with Keldysh 
Green's functions vanish in Eq. 

III. SUPERCURRENT 

We now investigate the current-phase relation /s((/)) in 
the parameter space which is controlled by temperature 
T and couphng strength T^^/A. This analysis focuses 
naturally on a calculation at finite T, for which regime 
the NCA yields quantitatively well controlled results^Z, in 
the absence of magnetic field, and so is expected to do 
its superconducting extension, the SNCA. We will com- 
pare the T ^ extrapolation of these calculations with 
our perturbative renormalization group (RG) analysis of 
the same model at T = as well as with exact numeri- 
cal renormalization group (NRG) calculations by Choi et 
al.'^'^ at T = for the symmetrical Anderson model. 

It has been elaborated by Clerk and Ambegaokar— 
that strong and weak coupling regimes are to be distin- 
guished by the position of the subgap resonance: the 
resonance moves through the Fermi energy from below 
when the coupling is increased, a behavior to be associ- 
ated with a transition from a tt- to a 0-junction type. In 
fact, we confirm this behavior in Figs. |21 and ^ How- 
ever, the resonance is wider than observed by Clerk and 
Ambegaokar: For strong coupling T/^/A > 1, the sub- 
gap resonance as well as the features at the gap edges in 
the d-electron spectrum are of the order of Tk > A, as 
Cooper pairs are broken in order to screen the impurity 
spin in this regime. The gap edges are less pronounced 
(see Fig. |31 lower panel), and a Fano-like interference 
between the continuum states and the subgap mode is 
evident. 



local spectral function 



150 




CO /A 

FIG. 3: The spectral function of the impurity d-level for weak 
coupling (upper panel, Tk/A — 0.125) and strong coupling 
(lower panel, Tk/A — 2.0). The leads are in the supercon- 
ducting state at T = 0.5A. The spectral functions are sen- 
sitive to the phase difference (j) = (j>L ~ (pa between left and 
right lead. 



The current-phase relation traverses three scenarios 
or transitions, as the coupling parameter Tk/A. is raised 
from weak to strong coupling (left column of Fig. 0}. 
These scenarios are related to the fact that the Josephson 
current states of a superconducting junction are equilib- 
rium states and are thus determined by the minima of 
the free energy. One may identify a succession of four 
current-carrying equilibrium states: 0-junction: single 
global minimum for (jj — 0; O'-junction: global minimum 
for (f) = and local minimum for (f> — n; 7r'-junction: local 
minimum for = and global minimum for (j) = tt; tt- 
junction: single global minimum for (j) = n. The succes- 
sion of the corresponding transitions has been discussed 
in the literature^2ii&. 
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FIG. 4: Phase dependence of the Josephson current for weak, 
intermediate and strong coupling values of Tk / A (left column 
of panels). In the right panel, the spectral function of the 
impurity d-level is presented for each coupling strength and 
temperature. 



1. Weak coupling. In the case Tk/^ — 1/8 all three 
curves /s(<^) correspond to temperatures above Tk- 
The curves are nearly identical within numerical 
resolution. We observe a 7r-junction behavior where 
the first harmonic, -sin0, dominates. In SNCA we 
cannot approach the low temperature limit where 
T is well below Tk- The latter has been investi- 
gated by Choi et al?"^ who indeed find a sinusoidal 
behavior for the zero temperature limit in the weak 
coupling regime. 

2. Intermediate coupling. Here, higher harmonics be- 
come important (as in the second row of Fig. 0| 
for Tk/A = 0.28). The derivative of the cur- 
rent at </) = changes sign with the temperature 
somewhere close to T = 0.4 A (middle left panel), 
which corresponds approximately to temperature 
and coupling where the spectral function crosses 
the Fermi energy. Choi et al.'^'' observe a discon- 
tinuous behavior in the current-phase relation for 
the intermediate coupling regime at zero tempera- 
ture however this discontinuity is smoothed for fi- 
nite temperature^. The distinction between the 
TT — tt', the tt' — 0' and the 0' — transitions, which 
classify the appearance and vanishing of the two 
minima of the free energy as mentioned above, is 
made by the characteristics of the current-phase 
relation: The sign reversal of the slope of Is{(f)) 
at vanishing cf) in the intermediate coupling regime 
signifies the tt' — tt transition. 

3. Strong coupling. The lower left panel in Fig. ^ 



shows the current-phase relation for A/Tk = 1. 
For this value of the coupling we are already in the 
strong coupling regime in the sense that the sub- 
gap resonance is clearly above the Fermi energy and 
the supercurrent is positive in the considered phase 
interval (0-junction behavior). The temperature of 
all curves is below Tk . The curve is more sinusoidal 
for the lowest temperature whereas it develops a 
flatter region for cj) close to tt for the higher tem- 
peratures. This is a precursor to the 0— tt transition 
and we will see below that the 7r-junction behavior 
may be recovered for higher temperature if Tk/A. 
is not too large. 

The phase diagram is now derived from the analysis of the 
extrema in the free energy: (i) the tt-tt' transition takes 
place when the maximum at phase (/> = turns into a 
local minimum of the free energy which is equivalent to 
the sign change of the slope of the current-phase relation 
for vanishing (j> (circles in Fig. [S)); (ii) the tt'-O' transi- 
tion refers to the point in the (T, rft:/A)-parameter space 
where the global minimum of the free energy switches 
from (p — TT to (f> = (triangles in Fig. O ; (iii) finally, 
the 0'~0 transition corresponds to the conversion of the 
local minimum at phase (f> = tt into a maximum, that is, 
the slope of Isi(t>) changes sign for phase tt (squares in 
Fig. El . 

The symmetrical BCS- Anderson model has also been 
investigated by Siano and Egger'*® using the quantum 
Monte Carlo technique. However, Choi et al. point 
out that Ref. does not consider the true low tem- 
perature limit and that the scales, such as the Kondo 
temperature, differ exponentially from the conventional 
definitions. 42*^ 

The zero temperature limit cannot be reached within 
the SNCA scheme. For this purpose we have performed 
a perturbative RG analysis, analogous to the poor man's 
scaling approach for the normal state^^. In the one- loop 
evaluation, the vertex from the impurity coupling term 



Jk 



Kondo 



(7) 



kk' 



generates particle-particle and particle-hole loops of con- 
duction and pseudofermion Green's functions. Here Jk 
is the Kondo coupling (J/f = y^/|e(j| for V oo) and 
t' (ct*) are the Pauli matrices in the impurity spin space 
(conduction electron spin space). 

While these diagrams renormalize the Kondo cou- 
pling in the normal state, one-loop contributions with 
an anomalous conduction electron propagator are to be 
included for the superconducting state. Although the 
corresponding vertex is missing in the bare Hamiltonian, 
the RG flow will generate the coupling which is of the 
form 



H - W^fV' ,f 



-ko 



(8) 
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FIG. 5: Phase diagram for the — tt transitions. As discussed 
in the text, the lower right area corresponds to the 0-phase 
and the upper left area to the 7r-phase. The SNCA data points 
refer to the tt — tt' transition (circles), the tt' — 0' transition 
(triangles), and the 0' — transition (squares). As seen from 
these curves, the respective transition points as a function of T 
scale roughly with logT^. The open inverse triangle presents 
the transition point from a spin doublet to a singlet state of 
a Kondo impurity in a bulk superconductor within the NRG 
evaluation of Satori et alm~. The star on the horizontal axis 
is approximately the transition point in the NRG analysis for 
the symmetric Anderson model of the quantum dot contact 
(Choi et alt~). The fat line at the horizontal axis is the 
regime where the perturbative RG analysis suggests the tt' — 0' 
transition at zero temperature. 



This coupling term will have the effect of cutting the RG 
flow for small T^/A, that is, in the perturbative regime. 
Only the coupling term with go2 will be renormalized 
under the RG flow, a consequence of spin conservation 
and the symmetry of the order parameter. 

With the initial conditions of isotropic spin coupling 
and zero potential scattering term, the following RG 
equations are obtainediS 



dJ 
d\nD 

dg 
dlnD 



= -Re 
Re 



D 



D 



- A2 



A 



3 A 



2-.g^ + --J 



D 



2D 



(9a) 
(9b) 



where D is half the band width and the dimensionless 
couplings are defined as: 



J — NqJk and g — ~iNo go2 



(10) 



Here g denotes the local coupling of pair fluctuations to 
the impurity with the initial condition g^Dg) = for 
the bare band cut-off Dq. For A ^ one recovers the 
standard poor man's scaling result. 



For weak coupling {Tk/A < 0.25) the effective band- 
width D approaches A before J or g diverge; the square 
root in the scaling Eqs. then vanishes and cuts off the 
RG flow. The scaling trajectories in the J~g plain flow 
towards a line of such fixed points. The ground state is 
an unscreened spin. The strong coupling regime is only 
accessible within a non-perturbative analysis but the ten- 
dency of the trajectories to flow away from the fixed point 
line towards a strong coupling fixed point with J — oo 
is already observed in the present one-loop evaluation. 
The strong coupling fixed point is approached for a bare 
coupling strength above a minimal Tx/A somewhere in 
between 0.25 and 0.65. 

The solid fat line at the horizontal axis in Fig. [S] in- 
dicates the range where this quantum phase transition 
(tt'-O' transition) is supposed to take place. The range 
is certainly too wide in order to estimate the low tem- 
perature extrapolation of the intermediate temperature 
data. However the transition range is consistent with the 
SCNA results. This does not apply for the extrapolated 
transition value of Clerk and Ambegaokar^^ which is al- 
ready at positive values of log(rf(-/A). We do not well 
understand the discrepancy to the result of Clerk and 
Ambegaokar; it may be related to the way in which the 
zero temperature limit was approached in Ref. [U How- 
ever the NRG result of Choi et al?'^ (with transition at 
Tk/A. ~ 0.42, star in Fig. O for the particle-hole sym- 
metric model is within our estimate from the 1-loop RG. 
Yet this comparison should be taken with caution as the 




0.2 0.4 0.6 0.8 1.0 

g 

FIG. 6: Scaling trajectories for various bare coupling values: 
(a) Tk/A = 0.1, (b) Tk/A = 0.2, (c) Tk/A = 0.5, (d) 
Tif/A = 1, and (e) T^/A = 2. The open circles represent 
fixed points, the dashed line is the fixed point line. The last 
fixed point in the perturbative regime is found for Ta /A ~ 
0.25. For Tif/A ~ 0.65 no downturn of the trajectory is 
observed. In the limit T^/A oo the trajectory follows the 
J-axis towards strong coupling. 
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particle-hole symmetric Anderson model also allows for 
local equal-time pair correlations (of the /-particles) and 
has no potential scattering term, both in contrast to the 
infinite-[/ case. The potential scattering term present 
in the asymmetric model induces a characteristic shift 
of the Kondo resonance relative to the Fermi energy^Si 
and thus is expected to influence the O-tt transition as 
well. The transition from the spin doublet to a singlet 
state of a Kondo impurity in a bulk superconductor was 
calculated within NRG by Satori et alm^. They found 
Tjf/A ~ 0.3 which is presented by the open triangle in 
Fig. El This NRG result is within the range of our 1-loop 
RG estimate and appears to be in agreement with the 
finite-temperature SNCA data. 

It should be obvious from the phase diagram that the 
transition is not only achieved by a change of the coupling 
parameter Tk/A (through, e.g., gating the quantum dot) 
but also by a temperature variation, provided that the 
(fixed) coupling is in an intermediate range. 
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IV. ZERO BIAS CONDUCTANCE 

The Josephson current directly probes the phase- 
sensitive anomalous Green's functions (see Eq. 0). The 
conductance, however, is related to the quasiparticle cur- 
rent and it is expressed through the imaginary parts of 
the diagonal Green's functions: it has a finite, mea- 
surable value if the derivative of the Fermi function 
is not exponentially small and if the bulk density of 
states (i.e., \mG^{Lj)) and the impurity spectral function 
0<^Gdi^)) both finite in the same frequency interval 
(see Eq. ©). This implies that the conductance van- 
ishes exponentially for the zero temperature limit as we 
restrict our considerations to s-wave superconductors. 

The question arises if this conductance of the quantum 
dot exhibits a signature of the O-tt transition at all. The 
solution to this elementary question is not straightfor- 
ward since the conductance does not directly expose the 
phase dependence of the superconducting states in the 
leads. However, as the position of the subgap resonance 
moves from below through the Fermi energy with increas- 
ing coupling, the ground state transits into a singlet state 
through the Kondo screening of the impurity spin. The 
enhanced screening, which is possible for strong coupling, 
not only modifies the subgap resonance but also the con- 
tinuum through increased pair breaking in this regime 
(cf. Fig. 13). Correspondingly, one may expect a feature 
in the temperature-dependent or coupling- (T/f/A-) de- 
pendent conductance which signifies the transition. This 
can be expected only if the temperature is not too small 
(with respect to A) so that df/duj in Eq. ^ is still siz- 
able for frequencies with finite p(w). 

Consequently, the conductance G, Eq. ©, represents 
an integral with respect to the quasiparticle energy lo 
over an interval of the order of T. In order to ex- 
hibit characteristic features of the local spectral density 
in the finite-T conductance it is, therefore, suggestive 



FIG. 7: Zero bias conductance G as a function of coupling 
strength Tk/A for fixed temperature. 



to analyse the temperature derivative dG/dT. Equiva- 
lently, one may consider the coupling parameter deriva- 
tive, dG/d{TK/A), because G is expected to be a uni- 
versal function in terms of T/Tk for fixed A. At finite T 
and fixed A, both dG/dT and dG/d{TK/A) may be ex- 
pected, in a rough first estimate, to be essentially propor- 
tional to \mQ^{uj = T) and \mQ^{u} = Tk), respectively, 
disregarding the energy dependence of the quasiparticle 
density of states in the leads — a more refined analy- 
ses certainly has to take the detailed frequency structure 
of the integrand into account. Hence, in the weak cou- 
pling regime (ln(A/rft:) > 1, ln(T/TK) > 1) we expect 



T — ' — I — ■ — I — ■ — I — ' — I — ' — r 




0.0 0.2 0.4 0.6 0.8 1.0 1.2 0.1 0.2 0.4 0.6 0.8 1 

T^/A 

FIG. 8: Derivative of the zero bias conductance, G{Tk / A), 
with respect to the coupling parameter at T/A = 0.36. Left 
frame: linear scales; the tt'-O' transition is at approximately 
Ta'/A ~ 0.4, for strong coupling dG/d{TK/A) is constant. 
Right frame: coupling scale is in logarithmic presentation; for 
weak coupling dG/d{TK/A) diverges logarithmically. 
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approximately that dG/d(TK/^) oc — ln(TR-/A), while 
in the strong coupling region (Tr-/A > 1, T/A < 1), 
dG/d{TK/A) should be approximately T and Tk inde- 
pendent. This expected linear behavior of G for large 
coupling strength, where Kondo screening is dominant, 
is associated with the formation of the Kondo resonance 
as it collects spectral weight and saturates at energies 
below Tfc- 

For the numerical evaluation we focus on the lowest 
temperature in Fig. [7| (T/A = 0.36, continuous line). 
The conductance displays a constant slope in the 0-phase 
{Tk/A > 0.4), whereas in the 7r-phase {Tr/A < 0.4) the 
slope of G tends to diverge for Tk 0. This crossover 
is even more apparent from the derivate of the conduc- 
tance with respect to Tk/A (see the left panel of Fig. |S| 
and the discussion above). Clear logarithmic behavior 
of the slope is observed in the 7r-phase, where the nu- 
merical SNCA results are especially controlled (see right 
panel of Fig. |SJ . The crossover between constant slope 
and logarithmic behavior is found to be at a coupling 
strength which correponds to the O-tt transition. In fact, 
for T/A = 0.36, the tt'-O' transition is at approximately 
Tx/A~ 0.4. 

Finally, we explore the temperature dependence of the 
quantum dot conductance. In Fig. we present the 
SNCA results for the zero bias conductance G versus 
temperature T at fixed values of Tk/A. The exponential 
regime for very low temperatures is outside the range 
where the SNCA is reliable. For the low temperature 
limit df /duj is exponentially small for w > | A| and should 
control the temperature dependence of the conductance. 
However, with our data we are still in the regime with 
a wide derivative of the Fermi function and a negative 
curvature of the G(T) Hues. Fig. El suggests to assign a 
constant slope to the conductance in this intermediate 



CD 




FIG. 9: Zero bias conductance G as a function of temperature 
T/A for fixed coupling strengths Tk/A. 



temperature range for the strong coupling regime. Such 
a linear behavior in G{T), Eq. is generated by an 
approximate compensation of the sharp spectral struc- 
tures in p{Lij) and in \mQ^{ijj) at the gap edge. The in- 
tegration over the derivative of the Fermi function yields 
the linear dependence for the intermediate temperature 
range — although the non-leading contributions to the 
frequency dependence above the gap edge may alter the 
temperature dependence. The linear temperature be- 
havior is consistent with the linear Tk dependence of 
G which was discussed before. For the weak coupling 
regime one should expect a InT dependence of dG/dT. 
The temperature range for the lowest curve in Fig. |2| 
(with Tk/A = 0.16) is too narrow to decide about a log- 
arithmic divergence of dG/dT but the few data points 
are consistent with this assumption. 



V. CONCLUSION 

In conclusion, we have calculated the Josephson cur- 
rent as well as the linear response quasiparticle conduc- 
tance for quantum dots in the Kondo regime with super- 
conducting leads. For finite temperatures, we have used 
a superconducting extension of the non-crossing approx- 
imation (SNCA), while for T = the behavior was an- 
alyzed using a perturbative renormalization group treat- 
ment. In this way we mapped out the phase diagram of 
the 0— TT transition of the Josephson current in the param- 
eter space of temperature and Kondo coupling constant 
for the first time in a systematical way. We stress that for 
temperatures not too far below Tk the SNCA is expected 
to produce reliable, semiquantitative results. The T — > 
extrapolation of our results agrees well with an NRG 
treatment of the problem at T = 0)2i considering that 
the latter was done for the &mte-U symmetric Anderson 
model which, on one hand, allows for equal-time Cooper 
pair formation on the dot, and, on the other hand, does 
not have a potential scattering term — both in contrast 
to our asymmetric, infinite-?/ model. Considering finite 
T will be essential for the analysis of experiments (see 
below) . 

Our results confirm that the Josephson current under- 
goes a succession of three transitions, — 0', 0' — tt', 
tt' — TT, separating four different Josephson equilibrium 
states, as the Kondo temperature Tk is reduced below 
the superconducting gap energy A, or, alternatively, as 
the temperature T is raised above Tk. In going from a 
0-type to a 7r-type junction, the four types of Josephson 
junctions are related to the successive development of the 
minima of the free energji^Si^. By explicit calculations 
we could relate these phases to the characteristic phase 
dependence of the Josephson current as well as to the po- 
sition of the Kondo- like subgap resonance above or below 
the Fermi energy. Moreover, we have also identified the 
signature of the O-tt transition in the quasiparticle linear 
response conductance G. Since the latter vanishes expo- 
nentially for T ^ 0, a treatment at finite T was essential 



here. From our results, the O-junction regime appears 
to be characterized by a constant slope of the conduc- 
tance as a function of Tr-ZA, while in the 7r-junction 
regime the slope diverges locarithmically. Note that at 
finite T the O-tt transitions are continuous crossovers be- 
cause of the finite width of the subgap resonance and 
develop a discontinuous jump only for T — > with van- 
ishing width of the subgap resonanceiSI These relations 
may be relevant for identifying and analyzing the dif- 
ferent phases in experiments like quantum dots, carbon 
nanotubes or other Kondo molecules coupled to super- 
conducting leads. Gated devices are supposed to control 
the level position in the quantum dot and, correspond- 
ingly, the Kondo scale. 
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APPENDIX A: SUPERCONDUCTING NCA 
(SNCA) 

In this appendix, we give a detailed derivation and dis- 
cussion of the formulae for the selfenergies, the auxiliary 
particle propagators and the local electron Green's func- 
tion within the SNCA. The same approximation has been 
used by Clerk and Ambegaokar in Refs. Ilj4l| 

The BCS- Anderson Hamiltonian Eq. obeys a local 
U(l) gauge symmetry with respect to simultaneous, time 
dependent transformations of the auxiliary particle fields. 



/^^e^^W/, , h 



which is intimately related to the conservation of the lo- 
cal charge Q, and which, due to Elitzur's theorem, cannot 
be broken. 54 The local gauge symmetry is implemented 
in the standard way^ where a conserving, selfconsistent 
approximation is generated from a Luttinger-Ward func- 
tional via functional derivative with respect to the renor- 
malized pseudoparticle propagators as well as the lead 
Green's functionsi^^*^ In addition, the constraint Q = 1 
is enforced in any expectation value of a physical oper- 
ator acting on the impurity state (more precisely: any 
operator which annihilates the |Q = 0) state) by taking 
the limit A oo; e.g. for the physical d-electron Green's 
functioni^Si^i 



Gdait) = -i lim 

A — >oo 



(d,(04(o) c-/^g-/^^(g-i) 

(g e-/3^^-/3A(Q-i)) 




FIG. 10; Pseudofermion selfenergy T.{iuj) and slave boson 
selfenergy II (ip) 



with (...) the time-ordered, grand canonical expectation 
value and /3 = 1/T. 

Implications of the projection Q — 1. The constraint 
crucially influences the auxiliary particle dynamics and, 
in particular, prohibits any anomalous contributions to 
the auxiliary particle propagators, even in the case of su- 
perconducting leads. Consider, for example, the Nambu 
pseudofermion propagator. 



(A3) 



(A2) 



c^-erf-A-E,(Lj) -S~(l^) 
-S— *(-^) c^ + ed + A + S_,(-c^) 

where Sg- is the normal selfenergy, and the anomalous 
selfenergy 5]°"°™ is assumed non-zero for the moment. 
Using the gauge (w — A) — *■ w, performing the matrix in- 
version in Eq. (|A3|) . and then taking the limit A — > cxd 
proves that all but the (11) element of the pseudofermion 
propagator, F^{ui) = [F^{uj)]ii, vanish. An analogous 
proof holds for the slave boson propagator. As a result, 
the (retarded) pseudofermion and slave boson propaga- 
tors F^{llj), B^iuf), respectively, have only normal con- 
tributions, 

F^ioj)^{oj-ed-X-Y:^{uj))~^ (A4a) 

B^ioj) = {oj-X- n^(t^))"^ , (A4b) 

where are the pseudofermion and slave 

boson selfenergies, respectively, and the spin index has 
been suppressed in the absence of a magnetic field. 

The evaluation of the Matsubara sum over a pseu- 
doparticle frequency shows that each closed pseudoparti- 
cle loop carries a fugacity factor e"'^'*'. Hence, upon the 
projection A ^ oo the Luttinger-Ward generating func- 
tional is comprised of diagrams which contain exactly one 
closed pseudoparticle loop: Its vanishing fugacity factor 
is cancelled by the corresponding factor e~^^ of the term 
{Q in the denominator of any physical 

expectation value [see Eq. (|A2|) ]. 

Definition of the SNCA . To define a conserving approx- 
imation that describes coherent Cooper pair transmission 
through a Kondo quantum dot, the generating functional 
shown diagrammatically in Fig. |21has been chosen. The 
first diagram of Fig. |21 represents the conventional non- 
crossing approximation and describes normal quasiparti- 
cle transmission. The second diagram of Fig. |21 contains 
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two crossing anomalous superconducting Green's func- 
tions in the leads. It is the simplest (i.e. lowest order in 
the hybridization Va) contribution to incorporate coher- 
ent Cooper pair tunneling in the Luttinger-Ward func- 
tional. Therefore, the corresponding approximation has 
been termed "superconducting non-crossing approxima- 
tion" (SNCA). The resulting auxiliary particle selfener- 
gies S, n, are shown in Fig. [TUland read, 

jde'f{e')gA^') 

aa' 

X B"{uj + e) B^iuj + e') F^{uj + e + e') , (A5a) 
n^H=E^ Jdempa{e)F%) 

a 

+ jd^fi^9a{e) jde'f(e')g:,{e') 

a' 

X F^-iuj + e) F'^iuj + e') B"^{uj + e + e') , (A5b) 

where /(e) is the Fermi distribution function and N = 2 
is the spin degeneracy. As seen from Fig. ^| the SNCA 
incorporates exactly two coherent Andreev transmission 
terms, one where a pseudofermion disappears from the 
dot and forms a Cooper pair in the superconducting 
lead, leaving an additional, virtual pseudofermion hole 
with opposite spin behind, and one describing the in- 
verse process. Incoherent, sequential Andreev processes 
are included in the propagators to infinite order via self- 
consistency. The set of equations (jA5|l can be further 
simplified in extracting the explicit phase-dependence of 
the conduction electron functions, i.e. 



propagators. This yields (see Fig.[TT)l. 

\mg^{Lu) = -J^ e^^" lmF^(e + lo) ImB^(e) , (A6a) 

Im^f (u.)= - r cos (0 J ^ e-P^j ^f{e') \gie')\ 

xlm [F^{e + uj)B^ie + e' + Lo)] 

xlm [F^{e + e') B^{e)] . (A6b) 

Without loss of generality we have set (t>L + <t>R — 0. Note 
that, although the local U(l) gauge symmetry on the dot 
prevents anomalous contributions to the auxiliary par- 
ticle propagators, anomalous physical electron Green's 
functions on the dot do exist. Physically this means that 
temporally retarded Cooper pairs on the quantum dot 
are indeed induced by the proximity effect, even though 
the formation of equal-time Cooper pairs is completely 
suppressed by the local Coulomb repulsion. 

Numerical evaluation of the SNCA. The Eqs. (|A4|) and 
(jA5|) form a closed set of non-linear integral equations for 
the auxiliary particle propagators which is solved numer- 
ically by iteration. The physical dot electron Green's 
functions, which determine the Josephson as well as the 
quasiparticle current (see Appendix B), are then com- 
puted using Eqs. (|A6|) . Note, however, that the Boltz- 
mann factors in Eqs. HA6|I strongly diverge for negative 
frequencies. Although this divergence is compensated by 
the threshold behavior of the pseudoparticle propagators, 
a numerical evaluation necessitates a re-formulation in 
terms of a new set of functions, lmF(a;), \mB{uj). We 
define them via the relations 

\mF^{u) = f{-uj)\mF{uj) , 
ImB^(w) = /(-w) \mB{uj) . 



a 

5]r,r,,5a(e)g:,(e') ^r^cos^ (^) m\W)\- 

aa' ^ 

Here, we have defined the phase difference (j) = (j)^ — 
and, and assumed symmetrical superconductors in 
the leads, |5rL(e)| = |ffi?(e)| = l5(e)l, Pde) = Pi?(e) = 
p{e) . By means of these relations it is seen immediately 
that only the terms with anomalous conduction electron 
propagators contribute to the phase dependence of the 
selfenergies in Eqs. ljA5|l . 

The equations for the normal and for the anomalous, 
physical dot electron Green's function are derived anal- 
ogously by functional derivative of the generating func- 
tional with respect to the corresponding lead electron 



Since the Boltzmann factors appear precisely in conjunc- 
tion with the integrals along the branch cuts of the auxil- 
iary particle Green's functions, it is possible to absorb all 
these exponentially diverging factors in \mF{Lu), \mB{uj) 
by observing e~^^ f{—Lu) = /(w). Details of this method 
as well as of the efficient treatment of the projection 
X oo can be found in Ref. i37t 




FIG. 11: Diagonal element Q^{iuj) and off-diagonal element 
J-^{iLu) of the impurity Green's function. 



11 



FIG. 12: Diagrams for G^^^ ^; The double line represents the 
dot electron Green's function, the single line the conduction 
electron green's function. The left diagram contributes to 
the normal current, the second diagram contributes to the 
supercurrent, Eq. lB4t 



APPENDIX B: JOSEPHSON CURRENT 
THROUGH AN INTERACTING REGION 

A formula for the supercurrent through an interacting 
regio n is de rived. We proceed along the line of refer- 
ences The charge current through the system can 
be expressed by means of the time derivative of the elec- 
tron numbers Ni^/^ in the left and right lead, i.e. 



Il/r - T{NL/R{t)) , 



(Bl) 



Note that the relation II — Ir holds. The time derivative 
is calculated using the Heisenberg equation of motion, 
yielding the expression 



(B2) 



with riL/}i = =f1- Using the definition for the lesser func- 
tion <^kCTa,CT'(^'^') = «(f^CT'(i')ck(T,a(i)), la may be rewrit- 
ten as 



2e v^,, ^ 



(B3) 



The diagrams contributing to G^^^ ^ are shown in figure 
[T^ In the following, the system is assumed to be in 



equilibrium, i.e. eV = 0. In this limit the normal current 
vanishes, and charge is transfered through the system 
only via the supercurrent Is, with 



2e 



"2 / — Re 

27r 



- Tl^"{u:)F^,{oj) . (B4) 



To derive the above equation we have used the relations 
G(w) = G<{u;) - G^H and G{u) = G«(w) - G<iuj) 
for the time ordered and anti-time ordered Green's func- 
tions, respectively, and G<{lu) = f{uj) (G^(w) - G^{uj)) 
for the lesser functions. For convenience we extract the 
explicit phase dependence from the off-diagonal Green's 
fimctions. 



where we have introduced the retarded, local conduction 
electron Green's function F^iuj) = X^k-^kal'^)- More- 
over, without loss of generality we have set 0i, -I- (pR — 0. 
Finally, using local charge conservation in the stationary 
case, Is^L = —Is,R = Is, obtain the following formula 
for the Josephson currenlii^ 



2e r 

Is = 7^ sin(0) / duj f{uj) 



h ttNq 



with r = Fl + Ft?. 



X Im 



(B5) 
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